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O; 

We derive amplitude equations for interface dynamics in pattern forming systems with long-range interactions. The 
basic condition for the applicability of the method developed here is that the bulk equations are linear and solvable 
by integral transforms. We arrive at the interface equation via long-wave asymptotics. As an example, we treat the 
Grinfeld instability, and we also give a result for the Saffman- Taylor instability. It turns out that the long-range 
interaction survives the long-wave limit and shows up in the final equation as a nonlocal and nonlinear term, a feature 
, that to our knowledge is not shared by any other known long-wave equation. The form of this particular equation 
£S) ' will then allow us to draw conclusions regarding the universal dynamics of systems in which nonlocal effects persist 
| at the level of the amplitude description. 
! PACS numbers: 05.10.-a, 04.25.-g, 45.70.Qj, 62.20.-x 
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O I. INTRODUCTION 

o 

A number of pattern formation problems involving an interface are computationally difficult despite the fact that 
the bulk equations are linear. Examples include Laplacian dynamics (diffusion-limited aggregation, limit of vanishing 
Peclet number in dendritic growth), flow problems in the Stokes approximation, and elastic problems. They are 
I ■ difficult because of the long-range self-interaction of the interface mediated by the external field. With the availability 
of modern computer power, these problems can usually be solved to good accuracy in two dimensions. But none of 
them has been successfully treated on a large scale in three dimensions so far. 

Whenever an interface equation can be derived, this offers the advantage of better tractability, both numerically and 
analytically. Interface equations for three-dimensional systems are two-dimensional, rendering much larger systems 
accessible than bulk simulations [EJ. For problems containing a length scale restricting the interaction distance, such 
as the diffusion length in solidification problems, interface equations are local in the long-wave limit [|, i.e., they 
are partial differential equations. The Laplace or Lame equations, however, do not have a cutoff distance. Hence, it is 
?3 | interesting and important to investigate how this feature affects the long-wave limit, if any. There have been interface 
equations containing nonlocal but linear terms in the literature before [Q, ||, As it will turn out, the nonlocal 
aspects of the system we consider are much more ferocious, leading to terms that are both nonlocal and nonlinear. 
We are also aware of the case of a previously derived amplitude equation containing nonlocal nonlinearities |7| ^|. 
However, the equation in question is weakly nonlinear whereas we will obtain a strongly nonlinear equation. 

Let us recall that close to the instability threshold for the merging of order, systems exhibiting a type I s instability 
in the nomenclature of Cross and Hohenberg |9j, i.e., an instability at finite wave number, are all expected to be 
described by a universal amplitude equation of the Ginzburg-Landau type. Otherwise, when the critical wave number 
approaches zero (type II S ) the dynamics is described by a long- wave equation such as the Kuramoto-Sivashinsky (KS) 
equation. A prerequisite for these derivations is that nonlocal interactions are screened off beyond some length (small 
in comparison with l/q c )- When long-range interactions (e.g. in electrostatic, magnetic, elastic systems, etc.) are 
effectively present, no universal equation has been known as yet for the latter case. The results of 0, |) pertain to the 
former case, i.e., they constitute a generalization of the Ginzburg-Landau description including nonlocal interactions. 
We derive here a universal form of dynamics in the long-wavelength limit when long-range interactions prevail and 
obtain a result that is complementary to the Kuramoto-Sivashinsky description. The derivation is exemplified on a 
system undergoing a surface instability due to elastic stress. Besides the universal feature of the resolved question, 
we present the practical virtue of this strategy. 
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Applying our approach to a particular system, a uniaxially strained solid undergoing the Grinfeld instability |T(], [TT| , 
a problem which has recently attracted the community of crystal growth [fl2| , [13[ |l4| , we obtain an equation that, 
when truncated to the leading-order nonlinearities should provide an example of the universal dynamics. As we shall 
see, this equation contains fewer parameters than to be expected which is most likely due to the fact that in the 
case of the Grinfeld instability the dynamics is variational. Introducing an additional free parameter, we are able to 
give the generic form of the universal equation. A derivation of this equation from symmetry considerations will be 
presented elsewhere p!5[ . 

The paper is organized as follows. In Sec. II, we write down the model equations for the elastically strained system 
and introduce the appropriate rescaling to perform an asymptotic analysis. Section III describes the asymptotic 
expansion and matching procedure and presents the final interface equation for two dynamical situations corresponding 
to conservative and to nonconservative dynamics. In Sec. IV, we give some simulation results for these equations. 
Section ^ contains a discussion of universality aspects and gives an example of a simulation of the generic equation. 



Finally, we will briefly summarize some conclusions and give an outlook in Sec. VI 



II. MODEL EQUATIONS AND NONDIMENSIONALIZATION 

The physics of the Grinfeld instability has been discussed in some detail in jl6| . We consider the case of a solid 
in contact with its melt; let us call the interface position £(x,t) (i.e., we restrict ourselves to two dimensions). It is 
related to the stress distribution via (lj, [l?], 18 : 



Ct(x,t) 



(1) 



Herein, 1/k is a mobility, p s the density of the solid, E is Young's modulus, v the Poisson number, 7 is the surface 

tension, Ap the density contrast between solid and liquid, g the gravitational acceleration, k — —Cxx/ (l + C^) 3 ^ 2 is 
the interface curvature, and partial derivatives of £ are denoted by subscripts, au and a nn are tangential and normal 
stresses at the interface, i.e., a nn = TiiCTijUj, ott = tiUijtj, where m and U are the components of the normal and 
the tangent to the interface. By <jq, we denote the externally imposed uniaxial stress, and the interface position is 
measured from the equilibrium position of a planar interface. (Setting k = and \a t t — &nn\ = \&o\ m (0) leads to 
a steady state £ = 0.) Equation (|l|) describes the mass current density which is proportional to chemical potential 
changes produced by strain, capillarity and gravity. 

To solve this equation, it must be supplemented with the equations describing mechanical equilibrium. Assuming 
linear isotropic elasticity, 

E 



Ukk 5ij , (2) 



lJ 1 + v V J l-2u 

where Uij is the strain tensor, expressible via the elastic displacements Ui according to Uij = 5 (diUj + djUi) , the 
condition for mechanical equilibrium, dj/Jij — (neglecting gravity in the bulk as a small cross effect), takes the form 
of the Lame equations: 

(l-2i^)V 2 u + V(V-u) = 0. (3) 

Note that the static version is sufficient since all interface motions are slow in comparison with sound propagation. 
Boundary conditions are 

cr„ n =cr„t = (4) 

at the interface, periodicity in the x direction with some prescribed wavelength 2n/q, and 

&ij — > oo SixSjx for z — > -00 . (5) 

Herein, a n t = ni<Jijtj, is the shear stress at the interface, which must evidently vanish if the latter is in contact 
with a fluid. Another cross effect has been neglected here, viz. that of capillary overpressure, leading to a curvature 
dependent jump of a nn across the interface. An overall constant pressure, not influencing the dynamics of the interface 
has been subtracted out of the definition of a nn . 

Due to the free boundary character, equations (|^) through (|5|) do not completely specify the problem and there is a 
need for an additional condition, which is just equation ([!]) expressing the fact that the chemical potential difference 
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between solid and liquid is proportional to the normal mass current across the interface. This means that we assume 
the interface to be microscopically rough (so that a linear relation holds between the current and its conjugate variable, 
the chemical potential difference) and the solid to be in contact with its liquid or vapor. 

Note however, that for the solution of the elastic problem at one instant of time with given interface position we do 
not need ([!]) (since we have neglected sound propagation effects). This suggests to attack the problem in two stages: 
first solve the elastic equations, then move the interface according to the elastic fields, and repeat the procedure for 
the new interface position. 

The question is how to obtain a closed interface equation, and to identify its generic form. In order to proceed, we 
first nondimensionalize the above equations by introducing two length scales l\ — jE/a 2 (l — v 2 ) (the Griffith length 
[19 1, apart from a factor of and 1% = (1 — ^ 2 )cr 2 /2A p gE (a gravity length) as well as a time scale r = kp 8 l\/^ 
[ IS | . Then we use Hooke's law (^) to express stresses by strains and reduce the latter by a dimensionless factor 
E~f(l + u)a . We find 



(u tt - u nn ) - 1 



+ K + a( > , 



(6) 



where a = I1/2I2, i-e. a > 0. A linear stability analysis of the full set of equations @ through (||) yields the dispersion 
relation f\M: 



w = 2\q\ 



(7) 



where u> and q are the nondimensional growth rate and wavenumber, respectively. The absolute value of q arises, 
because the eigensolutions of the elastic problem carry a factor exp(±iqx + \q\z), approaching zero for z — » — oo. A 
linear instability will arise for a < 1. 

As is obvious from the derivation, to be able to manifest itself the instability requires mass transport. Besides the 
case described of a liquid in contact with its melt, corresponding to nonconserved dynamics, another situation is of 
particular interest, due to possible applications in epitaxial growth. This is a solid in vacuum. Material transport will 
ordinarily be dominated by surface diffusion then. Consequently we have conserved dynamics, with eq. (^) replaced 
by 



CtOM) 



l 



d. 



1 



VT+Cl V^+Q V^+Cl 



1 r 
2 



(Utt 



1 



+ k + a£ 



(8) 



and the dispersion relation (Q) is simply multiplied by a factor q 2 on the right-hand side. Since the gravity term is 
normally negligible in experiments on epitaxial growth (unless a temperature gradient mimicks a strong gravitational 
field 0, |l7)), we may set a = in this case, obtaining a parameter-free equation. 

As it turns out, the asymptotic analysis can be performed together for both dynamical situations, since its most 
extensive part consists in solving the elastic problem (which is the same for both cases). Only afterwards is the 
solution inserted into the equation describing the interface evolution. At intermediate steps of the calculation, we 
will therefore write down equations for the nonconservative case only, but we shall give the final interface equation 
for both cases. 

Let us now consider a situation, where t\ is very large, i.e., we introduce a small parameter e, setting t\ = 0(1/ e). 
This corresponds to a very small stress, but we can nevertheless reach an unstable state by reducing gravity or the 
density difference, so ii becomes larger than l\j2. 

Then it will be natural to measure length scales in units of el\ [which is 0(1)] rather than l\ and time scales in units 
of e 2 T. Coordinates are transformed according to x = x/e, z = z/e, and t = t/e 2 . Referring to the new coordinates, 
the dispersion relation transforms into u> — 2e\q\ — q 2 — e 2 a. For small e, each term in this relation must behave as 
e 2 , hence we have for the wavenumber q ~ e. A long- wave equation should therefore be derivable. To perform the 
asymptotics, we introduce another set of coordinates via X — ex, Z = z — £(x,t), T — e 2 t, conveniently mapping 
the interface position to Z = 0. X and T are slow variables. We set U = u x — (1 — v)X , V = u z + v(Z + £), where 
u x and u z are the displacements, which makes the boundary conditions at — oo homogeneous (U — > 0, V — > 0). The 
transformed equations dl|,p|) read: 



<T(l + e 2 Ci) 3/2 = 
1 



1 - e 2 C|) (1 + X U - Cxd z U - dzV) + 2Cx [d z U + e 2 (d x V - Cxd z V)) 
\ (1 + 6 2 Ci) 2 - eQ xx {l + e 2 C 2 x ) 1/2 + eaC(l + ^xf) , 



(9) 
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(1 - 2u)d 2 z U + e 2 [2(1 - v) (d x - Q x dzf U + (d x - ( X 8 Z ) 3 z v\ = , (10) 

2(1 - v)d%V + (8 X - Cxdz) d z U + e 2 (l - 2v) (d x - t x d z f V = , (11) 

and the boundary conditions at the interface (Z = 0) become 

= (1 - v)d z V + vd x U - (1 - v)(xd z U 

+ e 2 [(1 - 2v) (C| - CxdxV) + (1 - (d Z V + d x U - (xd z ll)] , (12) 
= [d z U + e 2 (d x V - CxdzV)] (l - e 2 ( 2 x ) + 2e 2 ( x (d z V - d x U + Cxd z U - 1) . (13) 

Equations (^) through ( [l3|) constitute our starting point for the asymptotic analysis to be carried out. 

III. ASYMPTOTIC ANALYSIS 

These equations are singular in Z , i.e., the boundary conditions at infinity, where the implicit assumption Z = 0(1) 
does not hold anymore, can be satisfied only with the trivial solution C, x = 0. This difficulty usually does not appear 
with equations that have an internal length scale, in contrast to the elastic equations which are devoid of such a 
scale. To overcome it, one introduces an inner region, where Z = O(l), and an outer region, where rj = eZ = Oil) 
|^| , attempting then to match the solutions in their common domain of validity. In the inner region, Eqs. (|Io|) -(|l^) 
hold, and solving the bulk equations may be reduced to ordinary differential equations in Z via expansion in powers 
of e. In the outer region, derivatives with respect to r\ and X are the same order of magnitude. Therefore, true 
partial differential equations have to be solved. It is here that the linearity of the Lame equations becomes important, 
because it allows one to solve these equations in terms of Fourier transforms. Nontrivial aspects of the calculation 
then are largely due to the asymptotic matching of real space functions with Fourier transforms. 

Denoting field variables in the outer region by small letters, u(X, rj,T) = U (X, Z, T), v(X, n, T) — V(X, Z, T), we 
obtain the outer equations 

{l~2v)d 2 u + 2{l-v)(d x -eC x d, 1 ) 2 u + e(d x -eC X d v )d n v = 0, (14) 

(d x -eC X d n )d ri u + 2(l-is)ed 2 v + e{l-2v)(d x -e( x d ri ) 2 v = 0, (15) 



the inner ones being given by ( |10[|ll| ). Expanding the interface position and the fields in powers of e, we find that the 
expansion is singular in another sense. If it is assumed that C = 0(1), nonlinear terms will not arise in the interface 
equation, a fact that can also be inferred from symmetry considerations and power counting. Scalings that lead to a 
nontrivial result are £ = (9(e _1 ), u — 0(1), v = 0(e^ 1 ). We therefore write 

C(X, T) = e-^-i (X, T) + t (X, T) + <i (X,T) + ... (16) 

implying dz£(X,T) = 0(1). Next we set 

V(Z) = e- 1 V- 1 (Z) + V (Z) + eV 1 (Z) + ... = v(i 1 )=e- 1 v_ 1 (T 1 )+v (7 1 ) + ev 1 (T 1 ) + ..., 

(17) 

U(Z) = U (Z) + eU^Z) + . . . = u( V ) = u ( V ) + e Ul (v) + ■ ■ ■ , (18) 

where for brevity we have suppressed the X and T dependences. We then obtain as matching conditions (taking into 
account the e dependence of n) 

Uo(Z) -> u (0), (Z^oo), (19) 
V-x(Z) -> «_i(0), (Z -><»), (20) 
Ui(Z) - Zu' (0)+ui(0), (Z -*•(»), (21) 
V Q (Z) ~ Zv'_x(0) + v (0), (Z-foo), (22) 

where the prime denotes a derivative with respect to the argument (rj). To the two leading orders, the inner equations 
are solved by 

U a (X,Z,T) = B Q (X,T), (23) 
V^(X,Z,T) = D_ X (X,T), (24) 
U t (X,Z,T) = ZA 1 (X,T) + B 1 (X,T), (25) 
V Q (X,Z,T) = ZC (X,T) + D (X,T). (26) 
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The outer equations need to be solved only at lowest order in e. For brevity, we rename C-i(X, T) back to £(X,T). 
Confusion should not arise, since we will never need £o nor It should be kept in mind, however, that the new 
C(X, T) is 0(1), because it is the prefactor of e _1 . The outer solution then reads 

1 f°° 

u (X, V ,T) = — [aoiq^+boiq^ + aX^e^+^+^^dq, (27) 

7-oo 

1 f°° ( 3 - Av 
u_l(X,»j,T) = — / \-isign(q)[ao(q 7 T) + bo(q 7 T)(v + C(X,T))}+i b (q,T) 



2ir _ 

e i q X+\ q \( V+ C(X,T)) dq ^ (28) 

Using the boundary conditions (12,13]), we can express A\ and Co by Bq and D—i: 

= 7 {Cx [2(1 - v) + (1 - 2u)& + ((1 - v)& + 2 - v)B' Q {X)] 

(l-v)(l + Cx) 

- (l-i>-i>&)tf-i(*)} , (29) 

Co(^) = \ A2 -5 {Ci + - (BSW+CxW^'-iW)} , (30) 

(1 - i>) (1 + Cx) 

where we have left out the T dependence (which is "passively" present everywhere but important only in the final 
equation - the Lame equations are always solved at a fixed time). 

The matching conditions at leading order give four equations expressing Bq, Z?_i, A\, and Cp by Fourier integrals 
involving ao(q,T) and bo(q,T). For Bo and D—\, these integrals are essentially shown in Eqs. (|27j , |28| ) , we just need 
to set rj = there. The two remaining equations read: 

i too r -, 
A 1 (X) = — / a (q)\q\ + b (q)(l + \q\((X)) e^+^Wdq, (31) 

27T J -co 1 j 

Co(X) = {-igao(?)-i6o(?)[sigii(9)+gCW]+i8ign(q)(3-4i/)6o(g)} 

e^+^Wdg. (32) 

So we have six equations altogether, viz. Eqs. (|27j) through J32]), for the six unknown functions -Bo, -^l, Co, 

Oo, and &o [ Pol - Their solution determines the elastic fields, which then can be inserted into the interface equation 
(0). However, what we need is an explicit analytic solution. Its details are technical, so we just give some essential 
steps. First, we rewrite the equations in terms of Laplace transforms (rather: one-sided Fourier transforms) 

fio(a;) = i- ^ a (q)e i ^dq, (33) 
1 

2^ 

for example, 

B (X) = a (X -i((X)) + C(X)b (X -iCPO) + c.c. . (35) 

Setting w = X + i Z, 

2a (w) ee m(X,Z) +in{X,Z) , (36) 
26 (w) = o{X,Z) + ip(X,Z) , (37) 

with m, n, o, p real functions, we can then reduce our problem to just two (real) equations involving complex 
quantities. These are the Eqs. ( p9] , ^fj| ) with all quantities expressed by a' (w) and by bo(w), i.e., finally by to, n, o, 
and p. 

Introducing the abbreviations 

^)= ^-^ (1 - 2 f , (38) 

(i-^)(i + C x 2 ) 



b (x) EE -L / hy q :< ■<■<!,!. ( :!4! 



A 2 (X) 

X 5 (X) 
A 6 PO 



Cx(2^+(l-^)Cx 2 ) 

(i-.)(i + C* 2 ) 2 



1-1/-1/C 



(i -")(i + a 2 )' 

2 



c 



a-,) (i + Cx 2 Y 

-v+{l-v)tx 2 
(1-,) (l + Cx 2 ) 2 ' 
Cx {-v + (!-«/) Cx 2 ) 



(i- V )(i+cc a r 

these two equations become 

1 - X 2 {X)Cx + (3 - 4i/)A 3 (J0] o(X, -C) + [-MX) + X 3 (X)Cx] m(X, -Q 
= X^X) + (\ 2 (X) - \ 3 (X)Cx) (ox(X,-() 

(-1 + X 2 (X)Cx + As(X)) (n(X, -C) + Cpx(^, -C)) - 2(1 - 2v)X 5 {X) Cxp(X, -0 , 



-1 - A 5 (X) + A 6 (X)CxJ m(X, -0 + [-X 5 (X)Cx + (3 - 4z/)A 6 (X)J o(X, -C) 
= A 4 (X) - (-1 - A 5 (X) + X 6 (X)Cx) C o x (X, -0 + 2(1 - 2v) (l - X e (X)Cx) p(X, -0 
+ (x 5 (X)Cx + A 6 (X)) (n(X, -0 + (px(X, -0) , 



(39) 
(40) 
(41) 
(42) 
(43) 



(44) 



(45) 



where px and ox denote the derivatives of p(X, — 0, o(X, — with respect to X at fixed £. That is, we first take the 
derivative and only then allow £ = ((X). 

Since the four functions involved are real and imaginary parts of of analytic functions, only two of them are 
independent. Indeed, we have 



n(X, -0 

p(x, -0 



i 

7T _ 
7T 



(X'-X) m(X',0) 

(x' - xf + e 

{X' - X) o(X', 0) 



dX' 



dX' 



, (x'-xr+e 

Equations (Q) and (|45|), being linear, can be solved formally for o(X, — and m(X, — 0. This leads to 

o(X,-0 - r 1 (X)[n(X,-0 + CPx(X,-0] + r2(X)p(X,-0 + r 3 (X), 
m(X, -0 = r 4 (X) [n(X, -0 + C P x(*, -0] + r 5 (X)p(X, -0 + r 6 {X)Qo x (X, -0 
+ r 7 (X), 

where the functions r\ through 7-7 are given by 

n{x) 

r 2 (X) 



l-2v 



2C 



x 



(l-2i/)(l + &) 



r 3 (X) 
r 4 (X) 
r 5 (X) 



-c 



(l-2i/)(l + &) ' 

o, 

-1 + ^(1 + Cx 2 ) 



(46) 
(47) 

(48) 
(49) 

(50) 
(51) 

(52) 
(53) 
(54) 
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rn(X) 
r 7 (X) 



-1 , 



i + C 



(55) 
(56) 



x 



The computations leading from (]2^) through (32) to ( }48| , f49| ) are straightforward but heavy. Some intermediate steps 
(e.g., the simplification of the ri(X)f as well as a verification of the whole calculation at the end have been done using 
computer algebra (MAPLE). Due to the chosen scalings ( [l6| ) - (|l8|), Eqs. (48 4^) constitute an exact reformulation 
of the outer elastic problem, which is linear in u and v, so higher order equations have the same form as the leading 
order in e. To obtain analytic results, we have to employ some approximation. We assume that C(X), the expansion 
coefficient of e — 1 , is itself a numerically small quantity. Then expressions such as (^6 47) become Hilbert transforms. 
Denoting the Hilbert transform of f(X) by 



H 



(57) 



where the integral is to be taken as a principal value, as indicated by the bar, we find, as the limit ( < 1 of fTJ||) 

o(X,0) = -ri(X)n\m(X,0)] - r 2 {X)n\o(X,0)] +r 3 (X) , (58) 



m(X,0) 



r 5 (X)H o(X,0) +r 7 (X). 



(59) 



This set of linear equations can be solved iteratively. Assuming, and this is the second approximation, that ^x(X) is 
also numerically small, we may truncate the intcration after the first step, taking accurately into account only terms 

up to order ^(X). Using H H[f(X)] = — f(X), we can solve the truncated iteration to obtain explicit expressions 

for o(X, 0), m(X, 0) that read: 



o(X,0) = -(l-2is)r 3 (X) + H r r (X) 



m(X,0) = -(l-2v)r 7 {X)-2{l-v)(l-2v)H r s (X) 



(60) 
(61) 



Substituting o, m, n = — Tl[m], and p = — Ti.[o] back into the equations for A\, Co, Bq, and inserting these into 

(H) ~ (|H) an( i then back into we finally obtain the sought-for interface equation, in which we rename X into x 
and T into t (they are the same variables anyway): 



Ct(x,t) 

(i + C, 2 ) 1/2 




(i-a i + 



7T J [X 1 - X){1 + ( X { X ') 2 ) 



-dx' 



i + C, 2 



(i + C 2 ) 3/2 



(62) 



The arguments of £ (x or x' and t) have been suppressed wherever possible without ambiguity. 

The nonlocal term drives the instability. It is interesting to see whether ([32]) will reproduce a coarsening dynamics 
of the type discussed in [fl8|| . For in that case whether a groove grows or shrinks does not depend on local quantities 
(such as the curvature) alone. Hence, the scenario described cannot be expected from a merely local equation. 

This is the melt-crystallization case |i~8) . It is easily checked that linearization of this equation produces the linear 
dispersion relation (0). The corresponding equation with dynamics controlled by surface diffusion is obtained simply 

by preceding the right-hand side with the operator — [(1 + £ 2 ) -1 ^ 2 dx] '■ 



Ct(x,t) 



x x/i + Cf x \ 2 a + C, 2 ) 5 



(l-C 2 ) i + - -I — 

1 \ IT J (X'~X)(1 



Ux') 



+ Ux') 2 ) 



dx 1 



+ 



4C 2 ' 

i + C 2 . 



Cxx 

(i + C 2 ) 3/2 



(63) 



For the system described by this equation, coarsening dynamics does not yet seem to have been studied in the 
literature, therefore it will be briefly discussed in the next section. 
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IV. SOME SIMULATION RESULTS 



Both equations (|6^) and ( p3| ) have been simulated numerically. Figure [I] gives an example for period-doubling 
dynamics found in the nonconservative dynamics when an initially periodic sinusoidal structure is submitted to a 
small perturbation of twice its basic periodicity. A similar result is obtained for conservative dynamics and the same 




FIG. 1: Period doubling for the interface profile. The initial interface is a cosine with wavelength 2n and amplitude O.f , with 
every odd minimum made deeper by 5%. a — 0.97. Time interval between curves: At — 10.0. 



parameter a, see Fig. g. In both cases, a was close to 1, hence the instability was weak. For smaller values of a, we 
found stronger differences between the behavior of the two systems. The main statement that can be made, however, 
is that in both of them similar coarsening scenarios are operative, since coarsening allows a reduction of the overall 
elastic energy. Because of the restriction that the average interface height has to remain constant in the conservative 
case, the latter has coarsening proceed more slowly, as can be verified by comparison of Figs, [l] and ||. 




FIG. 2: The initial interface is a cosine with wavelength 2n and amplitude 0.1, with every odd minimum made deeper by 5%. 
a = 0.97. Time interval between curves: At = 10.0. 

In Fig. ||, we give the dynamics obtained from a local perturbation of a periodic interface. The process of coarsening 
is more complex than the simple period-doubling scenario suggested in |lq| . Frustration effects come into play when 
grooves compete via the nonlocal interaction, so it may happen that not the next-nearest neighbor of the leading 
groove is a winner of the competition but rather the third neighbor. 

An extensive overview of the properties of Eqs. (^,^3|) is not intended here. It will be the subject of a different 
publication pl| . Suffice it to say that we find these amplitude equations to reproduce several qualitative features of 
the dynamics of the system faithfully, including nonlocally induced coarsening. Differences with the full dynamics |l§[ ] 
appear at large amplitudes. For vanishing gravity (i.e., a = 0), the amplitude is not predicted to saturate in the full 
model but does so in the amplitude equation. This is not so surprising given that the latter was derived assuming 
small amplitudes (more precisely a small prefactor of the 0(e _1 ) term in the amplitude expansion). In principle, 
systematic improvement is possible by continuing the iteration scheme truncated after the first step and including 
higher powers of £ and Cx- In practice, this improvement results in additional technical difficulties. 
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FIG. 3: Coarsening of the interface profile. The fourth groove of the initial interface has been made deeper by 1%. a = 0.90. 
Initial wavelength 2n. Time interval between curves: At = 10.0. 

It is certainly not possible to claim genericity of Eqs. (^32||33]). Their derivation accounts for only the lowest-order 
(i.e., quadratic) nonlinearities exactly, while resumming some nonlinearities to infinite order, which are presumably 
specific for the elastic system. However, if there is a universal amplitude equation for the systems considered, then it 
must contain the terms that our derivation produces exactly, and it may be sufficient to consider just these terms. 

V. GENERIC EQUATIONS 

Before looking at the generic equation(s) derivable from our particular case, let us briefly recapitulate the general 
situation concerning amplitude equations of universal validity. 

There is an over-abundance of nonequilibrium systems which spontaneously build up an organized pattern from an 
initially structureless state when they are taken sufficiently far from equilibrium H . Typical examples are present in 
hydrodynamical systems (e.g., Rayleigh-Benard convection), chemical reactions (e.g., Turing systems), crystal growth 
(the Mullins-Sekerka instability), and so on. Despite the fact that the underlying physical and chemical mechanisms 
are diverse in these systems, sufficiently close to the instability threshold they all fall within the same universality 
class. Let q c denote the wavenumber of the emerging ordered pattern, and h(x,t) a field describing the pattern 
(which can stand for a component of the velocity field in hydrodynamics, or an interface position in crystal growth, 
etc..) along the x direction. Close to the threshold it is known [gj that h ~ A(x,t)e lq " x , where e lqcX describes the 
(rapid) periodic variation of the field due to the birth of order, while A is a slowly varying (slow with respect to e lqaX ) 
amplitude. A obeys the following equation 

A t = A + A xx - \A\ 2 A , (64) 

where coefficients can always be set to unity by an appropriate rescaling. This equation is universal in the sense 
that its form depends only on translational and rotational symmetries. It is usually referred to as the amplitude 
equation or the Ginzburg-Landau equation || . The amplitude equation describes the nature of the bifurcation from 
the structureless state as well as instabilities of the ordered pattern with respect to wavelength modulations (the 
Eckhaus instability 

A description in terms of a slowly varying amplitude makes a sense only if A varies in a slow fashion in comparison 
with the pattern wavelength l/q c . This requirement fails if q c — > 0. In that case a separation between a fast and slow 
variation is illegitimate. One has thus to resort to a (singular) expansion a la Sivashinsky Q]. Because q c is small, it 
provides an appropriate parameter of expansion. This is also known as the long-wavelength limit. In that limit the 
field h obeys generically the Kuramoto-Sivashinsky equation (or the damped form with A the damping rate), which 
is written in the canonical form, 

h t = -Xh - h xx - h xxxx + hi . (65) 

Again this equation is universal p2| , p3| in the limit q c — > 0. Because of the smallness of q c (long- wavelength) any 
nonlocal effect present in the original constitutive equations disappears, and one recovers again a local equation. 
This is true only because a length scale is present and serves as a cut-off. It is akin to the Debye length which 
makes the nonlocal response function of the Coulomb gas local. There is, however, a myriad of situations where 
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long-range interactions play a decisive role at all scales, and where it is hardly believable that the above equation (or 
a similar equation) should arise even in the long- wavelength limit. Typical situations with long range interactions are 
electrostatic and magnetic systems, problems with elastic fields, and so on. These systems are devoid of an intrinsic 
length-scale and it is highly desirable to derive the appropriate generic form of the corresponding evolution equation 
in the long- wavelength limit. A derivation on the basis of a grad ient expansion of an appropriately chosen general 
nonlinear operator in the spirit of p2f will be given elsewhere |15| . 

Here we argue that we can essentially read off the form of that equation from our result (|62|). To this end, we 
simply expand to second order in £ a and drop all higher-order terms. As noted before, the kept terms are exact to 
that order. The "generic" evolution equation with long range interaction obtained this way reads 

oo 

C* = -"C + Cxx-- / %^ds'-2£ 

7T J X' — X 

— oo 

Comparing this equation with the KS one (|6q), we see that it is not exactly a generalization. Since the nonlocal term 
drives the instability, the diffusive term need not be negative and hence there is no need for a fourth-order term. 
Rather, Eq. (|66| ) is complementary to the Kuramoto-Sivashinsky result. Of course, the analog of Eq. ( |6^ ) with a 
negative sign of C, xx may also appear and then the dispersion relation would need a different ultraviolet stabilization. 
If this term were fourth order in q, we would obtain a linear term oc C xxxx and an equation that is a common 
generalization of the KS equation and Eq. ( |66| ) . 

Closer inspection shows that this equation can be only a particular case of the generic equation. For in general the 
amplitude equation will have coefficients in front of each of its six terms. Four of these can be made equal to one via 
division by a common prefactor and rescaling of space, time and the amplitude £ itself. So our final equation should 
contain two independent nondimensional parameters, but it has only one, a. We suspect that this is due to a hidden 
symmetry, and a natural candidate for this symmetry is the variational nature of the basic evolution dynamics ([l]|]). 
In the fully generic equation, the additional parameter would be the ratio of the coefficients of the two nonlinear 
terms. It is therefore easy to conjecture the form of this equation: 

oo 

7T J X — X 

— oo 

where we have introduced a second parameter j3. The corresponding universal equation for conservative dynamics 
would be obtained by preceding the complete right-hand side of this equation by the operator —d x . 

We have simulated Eq. (^6|) and find that it still describes coarsening, this time without saturation. An example 
of a simulation is displayed in Fig. || We see a single groove or "finger" survive after a long time. 

On physical grounds, an non-saturating amplitude is to be expected in some situations involving the Grinfeld 
instability jl8). Nevertheless, even in such a case one should a priori not believe that a long- wave equation predicting 
an ever-increasing amplitude correctly describes the long-time asymptotics, because the interface eventually leaves 
the domain of validity of the long-wave assumption. This happens when the amplitude becomes larger than the 
wavelength. Therefore, an equation such as (|6?]) should be expected to describe some intermediate-time asymptotics 
at best. 

Things are different for other universal equations such as the Ginzburg-Landau equation (|^|) and the (damped) 
Kuramoto-Sivashinsky equation (|65|), which usually do work in the long-time limit. But this is due to the fact that in 
these cases the final (average) amplitude does saturate and is predicted to be small. However, if the cubic term in d&jj ) 
is positive, which necessitates a fifth-order describtion to obtain saturation of the amplitude, the final asymptotics 
cannot ordinarily be trusted either, because then the amplitude may not be small contrary to the original assumptions. 

On the other hand, while any amplitude description leading to non-saturating amplitudes may seem self-defeating 
in this sense, there appear to be systems, where the domain of validity of the amplitude equation extends beyond 
that of the assumptions made in its derivation. Recently, a singular long-wave equation was derived for step-flow 
growth without desorption [ p4[ . It predicts an amplitude growing as the square root of time, i.e., without bounds. 
Currently available evidence on the basis of numerical solutions of the exact model equations seems to corroborate 
that prediction. Hence, it is imperative to compare the long-time limit of the amplitude description to whatever 
information is available about the long-time behavior of the full system, in order to assess its validity. 

Equation ( |67j ) should be thought of as being the first in a hierarchy of generic equations describing the initial 
and intermediate stages of evolution of a number of different systems. It is relevant to systems the linear dispersion 
relation of which is quadratic. For other systems, such as the Saffman- Taylor system, where the dispersion relation 
is cubic, we will have a different linear term in the equation (the Laplacian of a Hilbert transform), and it is an 
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FIG. 4: Evolution according to the generic equation (|57|). a — 0.8, (3 = 1.0. The initial interface is a random structure (visible 
near z — 0). Time interval between curves: At = 1.0 up to t = 20.0, thereafter At = 0.1, and the evolution is shown up to 
t — 21.6. Note the spectacular acceleration of the "winning groove". 

interesting question whether the appearing nonlinearities will remain the same. Following the same strategy for the 
Saffman- Taylor problem, we arrive at the following equation 

oo oo 

Ct = _v r U^± dx , _i_B_f C^} dx , + VC 2 + lBUxxx t m 

7T J X — X I y X — X 

— oo — oo 

where V is the flat-interface velocity, 7 the surface tension, and B = 6 2 /12/i with b the thickness of the Hele Shaw cell 
and \x the fluid viscosity. Of course the linear terms are always quite trivial, but what is noteworthy is that f3 = 0, 
i.e., the nonlinear terms are local, like in the KS equation. (Under the usual assumptions for long-wave equations, 
the second local nonlinearity is small in comparison with the KS one.) We call this limit the Sivashinsky limit where 
nonlocality is present in the linear terms only, as was derived originally by Sivashinsky for flame propagation [Q. 
This contrasts with our equation (|67| ) which involves in the long-wavelength limit a nonlinear and nonlocal term. 
Therefore, this equation should introduce a new universality class with long range interactions. 



VI. CONCLUSIONS 



To summarize, our method of derivation is applicable whenever the general solution of the bulk equations can be 
found by a transformation method. Usually this means they must be linear (with constant coefficients). A second 
condition, requiring that the dispersion relation can be transformed to long-wave form, should be always satisfiable. 
The generalization to three dimensions is straightforward, in principle. Our approach opens up the road to a realm 
of amplitude equations constituting a hitherto unknown type, the investigation of which seems worthwhile both from 
the mathematical and physical points of view; a new line of research should be stimulated. 

The fact that the inner expansion produces only polynomials in Z (and hence the outer solution is identical to 
the uniform approximation p5|) strongly suggests that the derivation may even be achieved without the detour via 
a long- wave approach, using a regular perturbation expansion of the outer equations. This would seem to indicate 
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that we are not really dealing with an equation that is valid only in the long-wave limit but rather with an ordinary 
(but strongly nonlinear) amplitude equation valid at all wavelengths for small enough amplitudes, thus a posteriori 
justifying the formal device of rescaling to large wavelengths. On the other hand, the first of our approximations 
corresponds to setting \q\ £ equal to zero in the exponents of (|27| ) and (^8|), which, after taking into account all 
rescalings, means that the amplitude is small compared with the wavelength. 

Finally, the amplitude equation is much faster (and easier) to handle numerically than previous schemes. Like 
the KS equation, for example, its leading behavior described by Eq. ( |67[ ) should be generic for systems undergoing 
long- wavelength instabilities with long-range interactions. 
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